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^ ■ Abstract 

l> 

O 
l> . 

O ■ This paper describes a measurement of the Michel parameters, rj, p, ^, ^5, 

and the average Uj. helicity, h^y^, in r lepton decays together with the first 
measurement of the tensor coupling in the weak charged current. The t~^t~ 
pairs were produced at the LEP e^e~ collider at CERN from 1992 through 
1995 in the DELPHI detector. Assuming lepton universality in the decays of 



X 



D I the r the measured values of the parameters were: rj = —0.005 ± 0.036 ± 0.037, 

-fn ; p = 0.775 ± 0.023 ± 0.020, ^ = 0.929 ± 0.070 ± 0.030, ^6 = 0.779 ± 0.070 ± 0.028, 

hi,^ = — 0.997±0.027±0.011. The strength of the tensor coupling was measured 
K^'. to be K^ = —0.029 ±0.036 ±0.018. The first error is statistical and the second 

H I error is systematic in all cases. The results are consistent with the V — A 

structure of the weak charged current in decays of the r lepton. 



(Eur. Phys. J. C16(2000)229) 
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1 Introduction 

The Michel parameters [|I|, 77, p, ^, and C,6, are a set of experimentally accessible 
parameters which are bilinear combinations of ten complex coupling constants describing 
the couplings in the charged current decay of charged leptons. The Standard Model 
makes a specific prediction about the exact nature of the structure of the weak charged 
current, r leptons provide a unique environment in which to verify this prediction. Not 
only is the large mass of the r lepton (and thus an extensive range of decay channels) 
strong motivation to search for deviations from the Standard Model but the r also offers 
the possibility to test the hypothesis of lepton universality. 

The Michel parameters in r decays have been extensively studied by many experiments 
both at e+e" colliders running at the Z pole and at low energy machines 0,^ This paper 
describes an analysis of r decays using both the purely leptonic and the semi-leptonic 
(hadronic) decay modes, the latter being selected without any attempt to identify the 
specific decay channel. By grouping together all the semi-leptonic decays one can obtain a 
high efficiency and purity at the expense of a loss of sensitivity to the relevant parameters. 
This sensitivity is recuperated by splitting the semi-leptonic decay candidates into bins of 
invariant mass of the hadronic decay products, each bin being separately dominated by 
a different r decay mode. Results are presented both with and without the assumption 
of lepton universality. 

The measurement of the Michel parameters in the purely leptonic decay modes of the 
r allows limits to be placed on new physics. The large number of Michel parameters, how- 
ever, reduces the experimental sensitivity in placing these limits. Moreover, the Michel 
parameterisation does not cover the full variety of possible interactions; in particular it 
does not include terms with derivatives. However, a complementary test of a special 
type of new interaction is presented. In addition to testing new couplings with leptonic 
currents that conserve fermion chiralities, the possibility of an anomalous coupling of a 
leptonic charged tensor current is explored. 

2 The Michel parameters and i^r helicity 

The most general, lepton-number conserving, derivative free, local, Lorentz invariant 
four- lepton interaction matrix element, Ai, describing the leptonic decay r -^ IVi^r, 
(/ = e or yu), can be written as follows @J^J§]: 



V ^ N=S,V,T 



r 



N 



VViJn {{Vur)m\ ^N | Mr, 



which is characterised by spinors of definite chirality. G is a coupling constant, and the 
T N represent the various forms of the weak charged current allowed by Lorentz invariance. 
The n and m in Eqn. |l] are the chiralities of the neutrinos which are uniquely determined 
by a given iV, i and j. In the case of vector and axial- vector interactions the chirality 
of the neutrino is equal to the chirality of its associated charged lepton, while it is the 
opposite in the case of scalar, pseudoscalar and tensor interactions. In all cases we refer 
to the helicities and chiralities of particles; those of antiparticles are implicitly taken to 
have the opposite sign. 

The g^A parameters are complex coupling constants. There are 12 of these but, ex- 
cluding the possibility of the existence of a vector boson carrying a chiral charge, two of 



the constants, gj^^^ and gJifi, are identically zero. As the couplings can be complex, with 
an arbitrary phase, there are 19 independent parameters. The Standard Model V — A 
structure for the weak charged current predicts that gJ^L = 1 with all other couplings 
being identically zero. Neglecting phase space effects, the rate for the decay r~ -^ l^Ur^^i 
can be written [0,|| as 

, C^ml A 

^ ^ 1927r3 16' ^ ' 

with the definition 

A - <\9L\' + \9ln\' + \9k\' + \9L\')+n\9ln\' + \9lL\') 

+ n\9ln\' + \9U' + \91l\' + \9Il\') = 16- (3) 

From the above normalisation condition the maximum values that the coupling con- 
stants g^ can take are 2, 1 and l/\/3 for N = S,V and T respectively. The parameter 
G has to be measured from the decay rate and absorbs any deviation in the overall nor- 
malisation. The shapes of the spectra and the ratios of branching ratios, as used in this 
analysis, are insensitive to the overall normalisation and hence to G. 

The matrix element written in Eqn. [^ can be used to form the decay distribution of 
the leptonic r decay as follows: 

-— = Ho{xi)-VrH,{xi) (4) 

= [ho{xi) + r]hr,{xi) + php(xi)] - Vr[^h^{xi) + ^5h^5{xi)] (5) 

where xi = Ei/Emax is the normalised energy of the daughter lepton and Vr is the average 
r polarisation. Emax is the maximum kinematically allowed energy of the lepton, /. In 
the rest frame of the r, E^ax = 2m ' • ^^ ^^^ laboratory frame E^ax ~ Ej. or the beam 
energy Ef,(.am- The /I's at Born level are polynomials and are illustrated in Fig. |1|. The 
Michel parameters, rj, p, ^ and ^6, are bilinear combinations of the complex coupling 
constants |I]] and take the following form in terms of the complex coupling constants: 

V = lRe{6glj,gl*j, + 6g]^^g^^l + gf,j,gYl + gf,L9U + gln9U + gL9^*n)-, (6) 

P = m9lL? + M9U' + \9L? + \9L\' + \9k-29lL\' + \9lH-2glnn (7) 
C = -Wnn? + \9ln? - \9k\' - \9L?) + H\9ln? - \91l?) 

-i\9U - 3|^LP + 3|,LP - \9L\') + 4Re(,|,,Sl - ^lW^); (8) 

CS = m9lL\'-M9lH\' + \9lL\'-\9L\'+\9k-29lL\'-\9iR-29U'). (9) 

With the Standard Model predictions for these coupling constants the Michel parameters 
7], p, ^ and ^6 take on the values 0, |, 1 and | respectively. 

It is instructive to consider the physical significance of some of these parameters. A 
single measurement of p does not constrain the form of the interaction. For example, if 
p were to be measured to be |, as is the case for the Standard Model prediction, then 
this would not rule out any combination of the six couplings gf^^, gf^j^, gf^^, gf^^, g^p^, g\i^ 
with the other couplings being zero. Indeed a V+A structure would have a value of p of |. 
In this case one must examine the other parameters. For example, a V+A structure would 
mean that the parameter ^ would be equal to — 1 . The values of the Michel parameters 
for several examples of interaction types are given in Table |l]. 

The 7] parameter is of particular interest. It is sensitive to the low energy part of the 
decay lepton spectrum. It is practically impossible to measure r/ for r — > ez/^I^e decays 
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Figure 1: Polynomial functions in the laboratory frame for the r —>■ IviVr decay channel 
at Born level. The rj polynomial is normalised to 771^/771^-. 
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Table 1: The couplings and Michel parameter values for various mixtures of vector and 
axial- vector coupling at the two vertices in the decay r -^ Iv-rVi. 



because of a heavily suppressive factor of ^^ in the /z^ polynomial. This suppressive 
factor is of the order of ^ 1/17 for r -^ fJ^T^rT^n and hence all sensitivity to 7] is in this 
channel. The /i^ polynomial receives contributions from the interference between vector 
and scalar and vector aTid tcTisor interactions and is therefore particularly sensitive to 
non V — A interactions. If 77 7^ there would be two or more different couplings with 
opposite chiralities for the charged leptons and this would result in non-maximal parity 
and charge conjugation violation. In this case, ii V — A is assumed to be dominant, then 
the second coupling could be a Higgs type coupling with a right handed r and muon 0. 
The leptonic decay rates of the r lepton may be affected by the exchange of these 
non-standard charged scalar particles and these effects can be conveniently expressed 
through the parameter 7] W^F^. The generalised leptonic decay rate of the r becomes 
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where Gw is the coupling of the r to a lepton of type /, and equals the Fermi coupling 
constant if lepton universality holds. The functions / and g and the quantity rjjj^ are 



described in [|TT| . The parameter r'^^ is a factor due to electroweak radiative corrections, 
which to a good approximation has the value 0.9960 for both leptonic decay modes of the 

2 

r. The functions / and g are phase space factors. The factor /(— t) is equal to 1.0000 for 

2 

electrons and 0.9726 for muons. However, the function 4^^(7(^) equals 0.0012, whereas 

2 

the value of 4^f7(^) is relatively large, equal to 0.2168. Hence, under the assumption 
of lepton universality, a stringent limit on r/ in r — > /iz7^z/,- decays can be set on the basis 
of the branching ratio measurements, since to a good approximation (see discussion in 
section ITT) , 
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The variable P^ is defined as the probability that a right handed r will decay into a 
lepton of either handedness . This variable is related to the Michel parameters ^ and 
C,6 and to five of the complex coupling constants in the following way: 

= li^+'^^-m- (12) 

Hence the quantity P^ is a measure of the contributions of five coupling constants in- 
volving right handed r's. One can therefore see that measuring the parameters ^ and ^5 
is of considerable interest in studying the structure of the weak charged currents. 

The Michel parameters are restricted by boundary conditions. The leptonic decay 
rate of the r in Eqn. ^ has to be positive definite. Certain combinations of the Michel 
parameters lead to unphysical effects. It has been shown P^JTBIJT^ that the following 
constraints must be satisfied: 

< p < 1, (13) 

lel < 3, (14) 

p-m > 0, (15) 

9 -9p +17^5-3^1 > 0. (16) 

These inequalities describe the interior and surface of a tetrahedron in {p,^,^S) space. 
The first two conditions arise from the fact that the different couplings in the definitions 
of the Michel parameters occur in quadrature. The 3rd constraint can be found directly 
if the T decay rate in Eqn. ^ is forced to be positive definite for all values of xi. The 4th 
constraint is derived from the equations of two of the surfaces of the physically allowed 
tetrahedron. It is interesting to note that the Standard Model values of the Michel 
parameters are consistently at the edge of the allowed region (see Fig. ^. These relations 
are used in Section |^ to place limits on the coupling constants using the measured values 
of the Michel parameters. 

The decay width of the semi-leptonic decays of the r can be written, assuming vector 
and axial- vector couplings at the decay vertices as 

i£ = Ho{x)+VrH,{x) 

= Ho{x) - K^VrH,{x) (17) 

where a; is a polarisation sensitive variable in each decay channel. For the case of r — »• nur 
this variable is cos^*, the decay angle of the vr in the r rest frame, whilst for the two 



cases of r — i> pVj. and r ^>- aiv^- the variable used is the uo variable described in [|15]. The 
polarisation parameter h^^ is defined as 

2Re(t;^a;) 

V - -, f2~n 12' \^^! 

\Vr\ + |ar| 

where Vr and a^ are the vector and axial- vector couplings of the r lepton to W bosons. 
In the limit of a massless Vj. this is equivalent to the r neutrino helicity. 

Assuming that the boson exchanged in producing the t^t~ pair only involves vector 
and axial-vector type couplings then the helicities of the r"*" and r~ are almost 100% 
anti-correlated. This fact is used to construct the correlated spectra: 

^ "^'^ ^^'^'{H,{x,)-H,{x,)){H,{x,)-H,{x^)) + 



r dxidx2 2 

1-Vr 



(Hoix,) + H^ix^)) {H^ix^) + H^ix^)) (19) 



2 

in terms of the polarisation sensitive variable x (which is decay channel dependent). For 
leptonic decays the Hq and Hi functions are the polynomials described previously in 
Eqns. ^ and |l^ and the polarisation sensitive variable is the scaled energy xi. 

It follows therefore that by performing two-dimensional fits over this distribution one 
has full experimental access to all of the Michel parameters together with the r neutrino 
helicity, h^^, and the r polarisation, Vr-, with the caveat that only vector and axial- vector 
currents are assumed to contribute in the semileptonic decays. 

3 Anomalous tensor couplings 

The Lagrangian for the decay of the r can be written in the following way: 

where W^ is the weak charged current of the decay products of the W boson and k^ 
is a parameter which controls the strength of the tensor coupling. The choice of such 
a kind of interaction to test for the existence of new physics is inspired by experiments 
with semi-leptonic decays of pions ||16| and kaons [0], which show a deviation from the 
Standard Model which can be explained by the existence of an anomalous interaction with 
a tensor leptonic current |18|]. Since the new interaction explicitly contains derivatives. 



its effect on the distortion of the energy spectrum of charged leptons in r decays can not 
be described in terms of the known Michel parameters. Constraints will be placed on the 
parameter k^ from the analysis of both leptonic and semi-leptonic r decays, fixing the 
Michel parameters to their Standard Model values. The inclusion of the semi-leptonic 
channels significantly increases the sensitivity to the new tensor coupling and imposes 
stricter constraints. 

For purely leptonic decays, the matrix element takes the form 

M = -J={vi\^a\VVi)\{v^^\^a\Ur^)-i^^qf3{Vu^\aaf^\Urj,)j, (21) 

where q is the momentum of the W. The laboratory energy spectrum of the charged 
decay product can be expressed as 

— ^f{xi)+Vrg{xi), (22) 

dxi 



where xi is again the normahsed energy of the daughter lepton as defined in section ^. 
The expressions for f{xi) and g{xi), accounting for the new tensor interaction, were 



obtained in the rest frame of a decaying lepton |T^. Neglecting the mass of the final 
lepton and boosting along the r fiight direction gives 



f{xi) = 5-9xf + 4xf + 2«:f(l-xf), 
g{xi) = 1 - 9x2 + g^3 ^ 2k^{1 - 3xi + 2xf) 

These functions are shown in Fig. |[ 
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Figure 2: Polynomial functions for the tensor coupling contribution in the r —>■ IViUj- 
decay channel at Born level. Left plot shows f{xi) and right plot shows g{xi). In both 
plots the solid line illustrates the Standard Model case, k^ = 0, and the dashed line the 
case for k^ = 1. 



In the Born approximation the tensor interaction does not contribute to the process 
T -^ -KV^. Thus, among the main r semi-leptonic decay modes, only the decay modes 
T —>■ pUt- —>■ [2t:)vt and r — *■ aii/,- — ;► (Svr)//,- yield information about the new tensor 
interaction. The sensitivity can be increased by performing the analysis using the two 
angular variables 9* and ip, where 9* is the angle between the emitted final (pseudo) 
vector particle after a Lorentz transformation into the r rest frame and the inverse of 
the three-vector component of this Lorentz transformation. The angle if) is related to 
the angle of the p or ai decay products in the p or ai system and is sensitive to the 
polarisation of the hadronic system. These two variables are discussed in [jl5| . 



The r decay to a particle of spin 1 and mass nih and a neutrino has two amplitudes, 
Al and At-, representing longitudinal and transverse polarisation of the spin-1 particle 
respectively. From the expression for the decay helicity amplitudes. 
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where e* is the polarisation vector of the spin 1 particle with momentum q and helicity 
A, one obtains 

At _ V^rrih aT_ 

Al rrir ai ' 



(25) 



where a^ = 1 + i^^ 1'^ and ai = 1 + {m\/m^)K^ /2. Therefore for r — » (27r) u and 

r — !► (Svr) I/, 

d^iV 
oc (1 + Vr)H^ + (1 - P,)/f-, 



dcos6'* dcos'?/' 



where 



H^ = /lo(V') ('^L'^r COST] cos Y + ctT^^/i sin r/ sin Y 



+ /ii(^) 



[aLiTir smr/cos y 



if = /lo(V') (cti^^T COST] sin ^ — aTW/i sin r/ cos y 
+ hi{ip) f aim,- sin r/ sin Y + ctT^^/i cos ?7 cos Y 



Ot^/i cos ?7 sin y j + cty^^h sin^ y 
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Note that the angle rj used here is unrelated to to the Michel parameter of the same 
name. Neglecting terms of 0{m'^/E'^), the relation between 1] and 6* is 
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I — ml + (ml + ml) cos 6* 



m^ + m,1 + (?7i^ — m|) cos ^* 



(29) 



4 The DELPHI Detector 



The DELPHI detector is described in detail elsewhere [P0ipi| . The following is a 
summary of the subdetector units particularly relevant for this analysis. All these covered 
the full solid angle of the analysis except where specified. In the DELPHI reference frame 
the z-axis is taken along the direction of the e~ beam. The angle B is the polar angle 
defined with respect to the z-axis, (p is the azimuthal angle about this axis and r is the 
distance from this axis. The reconstruction of a charged particle trajectory in the barrel 
region of DELPHI resulted from a combination of the measurements in: 

• the Vertex Detector (VD), made of three layers of silicon micro-strip modules, at 
radii of 6.3, 9.0 and 11.0 cm from the beam axis. The space point precision in r-0 
was about 8 /xm, while the two track resolution was 100 fun. For the 1994 and 1995 
data the innermost and outermost layers of the VD were equipped with double sided 
silicon modules, giving two additional measurements of the z coordinate. 

• the Inner Detector (ID), with an inner radius of 12 cm and an outer radius of 28 
cm. A jet chamber measured 24 r-0 coordinates and provided track reconstruction. 
Its two track resolution in r-0 was 1 mm and its spatial precision 40 /im. It was 
surrounded by an outer part which served mainly for triggering purposes. This outer 
part was replaced for the 1995 data with a straw-tube detector containing much less 
material. 

• the Time Projection Chamber (TPC), extending from 30 cm to 122 cm in radius. 
This was the main detector for the track reconstruction. It provided up to 16 space 
points for pattern recognition and ionisation information extracted from 192 wires. 
Every 60° in there was a boundary region between read-out sectors about 1° wide 
which had no instrumentation. At cosG = there was a cathode plane which caused 



a reduced tracking efficiency in the polar angle range |cos0| < 0.035. The TPC had 
a two track resolution of about 1.5 cm in r-0 and in z. The measurement of the 
ionisation deposition had a typical precision of ±6%. 

• the Outer Detector (OD) with 5 layers of drift cells at a radius of 2 m from the beam 
axis, sandwiched between the RICH and HPC sub-detectors described below. Each 
layer provided a space point with 110 /im precision in r-0 and about 5 cm precision 
in z. 

These detectors were surrounded by a solenoidal magnet with a 1.2 Tesla field parallel 
to the z-axis. In addition to the detectors mentioned above, the identification of the r 
decay products relied on: 

• the barrel electromagnetic calorimeter, a High density Projection Chamber (HPC). 
This detector lay immediately outside the tracking detectors and inside the magnet 
coil. Eighteen radiation lengths deep for perpendicular incidence, its energy reso- 
lution was ^E/E = 0.31/E°-^4 © 0.027 where E is in units of GeV. It had a high 
granularity and provided a sampling of shower energies from nine layers in depth. It 
allowed a determination of the starting point of an electromagnetic shower with an 
accuracy of 0.6 mrad in polar angle and 3.1 mrad in azimuthal angle. The HPC had 
a modularity of 15° in azimuthal angle. Between modules there was a region with a 
width of about 1° in azimuth where the energy resolution was degraded. The HPC 
lay behind the OD and the Ring Imaging CHerenkov detector (RICH), not used in 
this analysis, which contained about 60% of a radiation length. 

• the Hadron CALorimeter (HCAL), sensitive to hadronic showers and minimum ion- 
ising particles. It was segmented in 4 layers in depth, with a granularity of 3.75° in 
polar angle and 2.96° in azimuthal angle. Lying outside the magnet solenoid, it had 
a depth of 110 cm of iron. 

• the barrel muon chambers consisting of two layers of drift chambers, the ffist one 
situated after 90 cm of iron and the second outside the hadron calorimeter. The 
acceptance in polar angle of the outer layer was slightly smaller than the other barrel 
detectors and covered the range |cosG| < 0.602. The polar angle range 0.602 < |cos0| 
was covered by the forward muon chambers in certain azimuthal zones. 

The DELPHI trigger was very efficient for r final states due to the redundancy existing 
between its different components. From the comparison of the response of independent 
components, a trigger efficiency of (99.98 ± 0.01)% has been derived. 

5 Particle identification and energy calibration 

The detector response was extensively studied using simulated data together with 
various test samples of real data where the identity of the particles was unambiguously 
known. Examples of such samples consisted of e"'"e~ -^ e"'"e~, and e"'"e~ — > n^fi~ events 
together with the radiative processes e^e~ -^ e^e~7 and e^e~ -^ fi^fi~'~f. Test samples 
using the redundancy of the detector were also used. An example of such a sample 
is r -^ 7r(n7r'^), (n > 0), selected by tagging the ii^ decay in the HPC. This sample 
was extensively used as a pure sample of charged hadrons to test the response of the 
calorimetry and muon chambers. 



5.1 TPC ionisation measurement 

The ionisation loss of a track as it travels through the TPC gives good separation 
between electrons and charged pions, particularly in the low momentum range. Because 
of the importance of this variable it was required that there were at least 28 anode wires 
used in the measurement. This reduced the sample by a small amount primarily due to 
particles being close to the boundary regions of the TPC sectors where a narrow non- 
instrumented strip was located. The dE/dx pull variable, Yi^^E/dxi ^°^ ^ particular particle 
hypothesis (j = e, yU, vr, K) is defined as 

^^'"'^/"^ ~ aidE/dx) ^ ^ 

where dE/dx^^^^ is the measured value, dE/dx ^{j) is the expected momentum depen- 
dent value for a hypothesis j and a{dE/dx) is the resolution of the measurement. 

5.2 Electromagnetic calorimetry 

The HPC is used for e, 7 and 7r° identification. For charged particles Eass is the energy 
deposited in the HPC. For electrons this energy should be (within experimental errors) 
equal to the measured value of the momentum. Muons, being minimum ionising particles, 
deposit only a small amount of energy in the calorimeter. Most charged hadrons interact 
deep in the HPC or in the HCAL and thus look like a minimum ionising particle in the 
early part or all of the HPC, with an increased energy deposition in the later layers if an 
interaction occurs in the HPC. 

The ratio of the energy deposition in the HPC to the reconstructed momentum has 
a peak at one for electrons and a rising distribution towards zero for hadrons. The pull 
variable, YIe/p^ is defined as 



where p' is the momentum refit without the use of the OD, described in Section 2^ below, 
and a{Eass/p'', Eass) is the expected resolution on Eass/p' for an electron with associated 
energy Eass- This variable gives particularly good separation at high momenta. 

5.3 Hadron calorimetry and muon identification 

The HCAL was used in particular for separating pions from muons. As muons travel 
through the HCAL they deposit a small amount of energy evenly through the 4 layers 
and travel on into the muon chambers whereas hadrons deposit all their energy late in the 
HPC and/or in the first layers of the HCAL so that they rarely penetrate through to the 
muon chambers. Therefore muons can be separated from hadrons by demanding energy 
associated to the particle in the last layer of the HCAL together with an associated hit 
in the muon chambers. To further distinguish muons from hadrons one can construct the 
variable E^iay, the average energy deposited in the HCAL per HCAL layer defined as: 

EMay = ^P^x sin' Q (32) 

-'■^layers 

where Ehcal is the total deposited energy in the HCAL; Niayers is the number of HCAL 
layers with an energy deposit and sin^ 6 smoothes out the angular dependence of the 
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energy response of the HCAL (see Fig. ^. This variable can be seen in Fig. 0. Note 
that the step behaviour around polar angles of 50° and 130° is due to the reduction in 
the number of layers hit in the HCAL where a niuon passes through a mixture of barrel 
geometry and end-cap geometry. 
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Figure 3: The HCAL response to muons (left plot) together with the variable Ehiay (right 
plot) for a sample of hadrons and muons in 1994 data (barrel region only). The crosses 
are the data, the solid histogram is the simulated sum of hadrons and muons and the 
hatched area is the simulated muons. 



5.4 Momentum determination and scale 

A good knowledge of the momentum and energy of charged particles is required for 
a Michel parameter analysis. This is especially true for the leptonic channels. As al- 
ready mentioned the momentum is measured by tracking the particles in a magnetic 
field as they traverse the detector. The precision on the component of momentum 
transverse to the beam direction, pt, obtained with the DELPHI tracking detectors was 
A{l/pt) = 0.0008(GeV/c)~^ for particles (except electrons) with the same momentum as 
the beam. Calibration of the momentum is performed with e"'"e~ -^ /^^/^~ events. For 
lower momenta the masses of the K^ and the J/ip are reconstructed to give an absolute 
momentum scale for particles other than electrons estimated, to a precision of 0.2% over 
the full momentum range. 

The determination of the momentum of electrons is more complicated. In passing 
through the RICH from the TPC to the OD, particles traverse about 60% of a radia- 
tion length. A large fraction of electrons therefore lose a substantial amount of energy 
through bremsstrahlung before they reach the OD. Due to this the standard momentum 
measurement of electrons would always tend to be biased to lower values. This effect is 
somewhat reduced through only using the measured momentum without using the OD, 
p'. The result is that this "refit momentum" shows a more Gaussian behaviour than the 
standard momentum fit. The best estimate for the momentum of the electron, p^i, is 
constructed in such a way as to benefit from the better resolution of the momentum mea- 
surement at low momentum and the smaller bremsstrahlung bias of the electromagnetic 
energy measurement. The reconstructed momentum and the electromagnetic energy were 
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combined through a weighted average which took into account the downward biases of 
the two respective measurements. The energy of the radiated photons was also added to 
the electromagnetic energy measurement to reduce further the effects of bremsstrahlung. 

An algorithm was used which performed a weighted average depending on the value of 
Eass/p' ■ The further this value was from unity, the more the weight of the estimator with 
the lower value was down scaled relative to the other. The scaling factor was inversely 
proportional to the square of the number of standard deviations by which the value of 
Easslv' differed from unity. 

Subsequent references to the momenta of electrons imply the use of the best estimator 
Pel- The momenta of other particles are measured using the standard momentum fit, p, 
of the particle as it traverses the detector. 

6 The selection of the event sample 

In order to determine the Michel parameters, a sample of exclusively selected leptonic 
decays of the r together with an inclusive sample of semi-leptonic decays have been used. 
The data sample corresponds to the data taken by DELPHI during 1992 (22.9 ph~^ at 
Eem = 91.3 GeV), 1993 (15.7 pb'^ at E^m = 91.2 GeV, 9.4 pb-^ at E^n, = 89.2 GeV and 
4.5 pb-^ at Eem = 93.2 GeV), 1994 (47.4 pb'^ at E^m = 91.2 GeV) and 1995 (14.3 pb~^ 
at Ecm = 91.2 GeV, 9.2 pb~^ at E^m = 89.2 GeV and 9.3 pb'^ at E^m = 93.2 GeV). 

In all analyses, samples of simulated events were used which had been passed through 
a detailed simulation of the detector response |^ and reconstructed with the same pro- 



gram as the real data. The Monte Carlo event generators used were: KORALZ 4.0 p2|] 



together with the TAUOLA 2.5 EH] r decay package for e"'"e -^ t^t events; 



DYMU3 [g^ for e+e" -^ /x+zx" events; BABAMC [g|] for e+e" -^ e+e" events; JETSET 



7.3 [^ for e+e — > qq events; Berends-Daverveldt-Kleiss |2^ for e"'"e — >■ e+e e+e 
e+e~ — > e+e~/i+/i^ and e+e~ — > e'^e~T'^T~ events; TWOGAM |^ for e+e~ — > e+e~qq 



events. 

The variables used in the initial preselection of the r sample together with the selection 
of the various decay channels are described below. 

6.1 The e+e" — > t'^t~ sample 

At LEP energies, a t^t~ event appears as two highly collimated low multiplicity jets 
in approximately opposite directions. An event was separated into hemispheres by a 
plane perpendicular to the event thrust axis, where the thrust was calculated using all 
charged particles with momentum greater than 0.6 GeV/c. To be included in the sample, 
it was required that the highest momentum charged particle in at least one of the two 
hemispheres lie in the polar angle range | cosG| < 0.732. 

Background from e+e~ — > qq events was reduced by requiring a charged particle mul- 
tiplicity less than six and a minimum thrust value of 0.996. The e+e~ -^ qq background 
is however negligible in the analysis of the Michel parameters as one is looking for events 
with only one charged particle in each hemisphere. 

Cosmic rays and beam gas interactions were rejected by requiring that the highest 
momentum charged particle in each hemisphere have a point of closest approach to the 
interaction region less than 4.5 cm in z and less than 1.5 cm in the r — plane. It 
was furthermore required that these particles have a difference in z of their points of 
closest approach at the interaction region of less than 3 cm. The offset in z of tracks 
in opposite hemispheres of the TPC was sensitive to the time of passage of a cosmic 
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ray event with respect to the interaction time of the beams. The background left in the 
selected sample was computed from the data by interpolating the distributions outside 
the selected regions. 

Two-photon events were removed by requiring a total energy in the event, Ems, greater 
than 8 GeV and a total transverse component of the vector sum of the charged particle 
momenta in the event, p™*'"^, greater than 0.4 GeV/c. 

Contamination from e^e~ -^ e^e~ and e"'"e~ -^ ^^fi~ events was reduced by requiring 
that the event acoUinearity, 6acoi = cos~'^(— i^^'Pi ), be greater than 0.5°. The variables pi 
and p2 are the momenta of the highest momenta charged particles in hemisphere 1 and 
2 respectively. 

The e"'"e~ — >• e"'"e" background is reduced in the second instance with a cut on the 



radial energy Erad (defined as Erad = y Ei"^ + E2^ / Ebeam where Ei and E2 are the energies 
deposited in the HPC in a 30° cone around the highest momentum charged particle in 
each hemisphere and Ebeam is the beam energy). Events are retained if Erad < 1- 

The e"'"e~ -^ /i"*"//" background is reduced in the second instance with a cut on the 
radial momentum p^ad (defined as Prad = \/pi^ + P2^ /Pbeam where pi and p2 are the 
momenta of the highest momentum charged particles in each hemisphere and Pbeam is the 
beam momentum). Cutting on this quantity is also effective in reducing the e"'"e~ — >■ e"'"e~ 
background. Events are retained ii Prad < 1- 

As a result of the above selection ~ 93000 e'^e^ -^ t^t^ candidates were selected from 
the 1992 to 1995 data set. The efficiency of selection in the 47r solid angle was ~ 54%. The 
background arising from e+e" -^ e+e" events was estimated to be (1.07 ± 0.32)%, from 
e"'"e~ — > yU^/i" events (0.30 ±0.09)% and from four-fermion processes (0.93 ±0.28)%. The 
background from e^e~ — *> e^e~qq was negligible. Since the efficiencies and backgrounds 
varied slightly from year to year the data sets were treated independently. 

6.2 The T — ^ eFgi^r channel 

The r — »• ez/e^'r decay has the signature of an isolated charged particle which pro- 
duces an electromagnetic shower in the calorimetry. The produced electrons are ultra- 
relativistic and leave an ionisation deposition in the Time Projection Chamber corre- 
sponding to the plateau region above the relativistic rise. Backgrounds from other r 
decays arise principally from one-prong hadronic decays where either the hadron inter- 
acts early in the electromagnetic calorimeter or an accompanying 7r° decay is wrongly 
associated to the charged particle track. 

As an initial step in electron identification it was required that there be one charged 
particle in the hemisphere with a momentum greater than O.Olpfeeam- To ensure optimal 
use of the HPC it was required that the track lie in the polar angle range 0.035 < |cos G| < 
0.707 and that the track extrapolation to the HPC should lie outside any HPC azimuthal 
boundary region, as described in Section Ij. 

The dE/dx measurement is crucial to the analysis and so it was required that there 
were at least 28 anode wires with ionisation information in the TPC. It was required 
that the dE/dx measurement be consistent with that of an electron by requiring that 
the YY" dE/dx variable be greater than -2. This requirement was efficient, especially at low 
momentum, in retaining signal and removing backgrounds from muons and hadrons. 

The selection continued with a logical "OR" of two criteria, the first on the YV dE/dx 
variable, which was particularly good at low momentum, and the second on the YIe/p 
variable, which was particularly good at high momentum. A particle was taken to be an 
electron if it deposited greater than 0.5 GeV in the HPC and the value of YIe/p was greater 
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than -2 "OR" the measured value of YVdE/dx was greater than 3 and the momentum was 
greater than O.Olpbeam- The "OR" thus gives a high constant efficiency over the whole 
momentum range. The YV'dE/dx ci-nd YIe/p variables can be seen in Fig. ^. 

DELPHI 




dE/dx 



Figure 4: The IT dE/dx ^^^ Ue/p variables after application of all the other selection cuts 
except the one shown for 1994 data. The crosses are the data, the solid histogram is the 
sum of the signal and background and the shaded area is the background from r -/->■ evV 
events. 

The remaining background was reduced by requiring that there be no hits in the muon 
chambers and no deposited energy beyond the first HCAL layer. Residual background 



from r 



TTimr 



V was reduced by cutting on the energy of the most energetic neutral 



shower in the HPC observed in an 18° cone around the track. Neutral showers were not 
included in this requirement if they were within 1° of the track and hence compatible 
with being bremsstrahlung photons. 

The identification criteria were studied on test samples of real data. The efficiency of 
the dE/dx and HPC cuts were tested across the whole momentum range by exploiting 
the redundancy of the two. Since the simulation showed that the two measurements were 
instrumentally uncorrelated, the overall bin by bin efficiency was calculated from these 
two independent measurements. 

Backgrounds arising from non-r sources consisted of e^e^ — >■ e^e^ and four-fermion 
e"'"e~ — >• e"'"e~e"'"e~ events. The e"'"e~ -^ e"'"e~ background was suppressed by the standard 
r preselection cuts, i.e. Prad < 1 and E^ad < 1- Four-fermion events remaining after the 
E^is and P™*'^ cuts were further suppressed by demanding that if the r -^ ez/T-z/g candidate 
had a momentum less than 0.2i?{,eam and there was only one particle detected in the 
opposite hemisphere with similarly a momentum below 0.2i?fteam then the r — >■ euT-u^. 
candidate was retained if YYdE/dx fo'^ the particle in the opposite hemisphere was less than 
3 and therefore inconsistent with being an electron. 

Application of the above procedure on the 1992 to 1995 data resulted in a sample 
of ~ 21500 r -^ eVeVT- candidates. The efficiency of selection within the 47r angular 
acceptance was 35%. The background arising from r -/^ eV^v-r processes was estimated to 
be (3.89 ± 1.17)%, from e"'"e~ — >• e"'"e~ events (1.61 ± 0.48)% and from e"'"e~ — >• e"'"e~e"'"e~ 
events (0.53 ± 0.16)%. 



14 



6.3 The T — ^ /xi/^1/7- channel 

A muon candidate in the decay r -^ hT'^i't appears as a minimum ionising particle in 
the hadron calorimeter, penetrating through to the muon chambers. Due to ionisation 
loss, a minimum momentum of about 2 GeV/c is required for a muon to pass through 
the hadron calorimeter and into the muon chambers. It was therefore required that there 
be one charged particle in the hemisphere with sufficient energy to penetrate through the 
detector into the muon chambers. The candidate had to have a momentum greater than 
0.05pbeam and lie within the polar angle interval 0.035 < |cos9| < 0.732. 

Positive muon identification required that the particle deposited energy deep in the 
HCAL or had a hit in the muon chambers. This was achieved specifically in the first 
instance by insisting that the average energy per HCAL layer Ehiay be less than 2 GeV. 
A logical "OR" of two variables was also used in the selection. The track was required to 
either have a maximum deposited energy in any HCAL layer of less than 3 GeV together 
with deposited energy greater than 0.2 GeV in the last HCAL layer, or have at least one 
hit in the muon chambers. This combination of cuts gave a reasonably constant efficiency 
over the whole momentum range. The two selection variables, the energy deposited in 
the last HCAL layer and the number of hits in the muon chambers, can be seen in Fig. ^. 
The background was suppressed further by requiring that the sum of the energies of 
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Figure 5: The number of hits in the muon chambers and the energy deposited in the last 
layer of the HCAL after application of all the other selection cuts except the one shown 
for 1993 data. The crosses are the data, the solid histogram is the sum in simulation 
of the signal and background and the shaded area is the background from r -/->■ ^Vr^^ 
events. 



all the electromagnetic neutral showers in an 18° cone around the track did not exceed 
2 GeV. This cut was effective in further suppressing r -^ 7T{mr^) and e^e~ —>■ fi'^fi~'~f 
events. 

The identification criteria were studied on test samples of real data. The efficiencies 
of the HCAL and muon chamber cuts were tested across the whole momentum range 
by exploiting the redundancy of the two. After correcting the simulated data for a 
discrepancy in the depth of the energy deposition by hadrons in the HCAL the data were 
found to be well described. 
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Backgrounds arising from non-r sources consisted mainly of e^e^ —* f^^fJ'~, 
e"'"e~ — *> e"'"e~/i+/i~, e+e~ — >• e"'"e~r"'"r~ and cosmic ray events. The e"'"e~ — >• /I'^/i" back- 
ground was suppressed by the standard preselection cut, i.e. prad < 1- The remaining 
background was further suppressed by demanding that the event was rejected if there 
was an identified muon in each hemisphere with momentum greater than O.SEbeam and 
the total visible energy was greater than 70% of the centre-of-mass energy. The event 
was also rejected if the momentum of the identified muon was greater than O.Spbeam and 
the momentum of the leading track in the opposite hemisphere was greater than O.Spbeam- 

The four-fermion events e"'"e~ — *> e~^e~^~^fi~ and e"'"e~ — >■ e"'"e^r+r^, although back- 
ground processes, required no further suppression. 

Candidate r —>■ /iz/^z/,- decays with muon momenta below 2 GeV/c were selected 
with different criteria. At these energies muons do not have sufficient energy to penetrate 
through the HCAL to reach the muon chambers, thus making the selection more difficult. 
Instead, at these lower momenta, muon candidates were selected if the particle was seen 
in the last 3 layers of the HCAL. This procedure was tested using a sample of hadrons 
selected from the data and simulation by tagging p decays through the presence of a 7r° in 
the HPC. In order to study the signal, various variables were compared in the data and 
simulated data to see if the simulation correctly modelled the performance of DELPHI at 
these low energies. The response of the HCAL to these hadrons and muons with momenta 
below 2 GeV/c was well described by the simulation, after the correction described above 
for the hadronic showers. 

As a result of the above procedure ~ 26000 r — > jmu^i't candidates were selected 
from the 1992 to 1995 data. The efficiency of selection within the 47r angular acceptance 
was 45%, the background arising from r -/^ [u/^Ur processes was estimated to be (1.88 ± 
0.56)%, from e"'"e~ — > /i+/i^ events (0.52±0.16)%, from e"'"e~ -^ e'^e~fi~^fi~ events (0.58± 
0.17)%, from e"'"e~ — >■ e"'"e~r"'"r~ events (0.48 ± 0.14)% and from cosmic-rays (0.14 ± 
0.04)%. 

6.4 The T — ^ h{mT^)vT channel 

The r -^ inclusive one-prong hadrons channel makes no distinction between the pri- 
mary semi-leptonic decays namely r -^ ttUt, t -^ pi^r and r — >■ aii^r- Instead each decay 
candidate is separated into bins of invariant mass, constructed from the 4-momenta of 
the charged particles and all reconstructed photons. The invariant mass bins used were 
Miny < 0.3 GeV/c2, 0.3 GeV/c^ < Mi„^ < 0.95 GeV/c^ and M^^^ > 0.95 GeV/c^. 

The preselection of the r's for this channel is slightly different to that for the lep- 
tonic channels due to the smaller potential backgrounds arising from di-lepton events. 
Therefore there is no Prad cut in the preselection and the E^ad cut is loosened to 1.1. 

In order to identify hadrons one is forced to use almost all the components of the 
detector. To be identified as a hadron it was required that one particle was detected in 
a given hemisphere in the angular range 0.035 < |cosB| < 0.732. In the case of more 
than one particle being detected, the hemisphere was retained if the highest momentum 
particle was the only particle having associated vertex detector hits. This ensured that 
one also retains a high efficiency for one-prong r decays containing conversions within 
the detector. 

Further cuts were made depending on the invariant mass of the decay products. Fig. |^ 
shows the invariant mass distribution for all preselected r's, calculated assuming that 
all charged particles were pions and all neutrals were photons. Most background from 
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leptons comes at low invariant mass. Hence one should apply stricter criteria for these 
events. 

DELPHI 




GeV/c' 



Figure 6: The invariant mass distribution for all preselected r decays in 1995 data. 
The crosses are the real data, the solid histogram is the simulated e"'"e~ -^ t^t~ data 
together with the simulated background, the shaded area is the sum of the e"'"e~ -^ e^e~ , 
e"'"e~ — * /x^/x~ and the leptonic r decays. The pole at the n"^ mass is not plotted. 



The background from electrons was suppressed with the following two cuts. Firstly 
the measured dE/dx in the TPC had to be consistent with being a pion, so YV dE/dx < 2. 
Because of the importance of the dE/dx measurement to the selection it was also required 
that there were at least 28 anode wires with an ionisation measurement. This cut is 
particularly effective at low momentum. 

The second cut required that either the particle deposited an energy beyond the first 
layer of the HCAL or that the associated energy in the first four layers of the HPC be 
less than 1 GeV for invariant masses below 0.3 GeV/c^, and 5 GeV otherwise. This cut 
is particularly effective at high momentum. The combination of the two cuts therefore 
leads to an even efficiency for the suppression of electrons across the whole momentum 
range. 

Rejection of background from muons was only performed for events with invariant 
masses less than 0.3 GeV/c^. Muon background in higher invariant mass bins was found 
to be small enough to justify no further suppression. The muon rejection was based on 
the average energy per HCAL layer, E^iay It was required that either E^iay was greater 
than 2 GeV or that there was no energy deposited in the HCAL. In addition to this 
criterion it was also required that there were no hits in the muon chambers and that the 
momentum of the leading charged particle was greater than O.Obpbeam in order that it 
had sufficient energy to reach the muon chambers. For regions not covered by the muon 
chambers it was required that there was no deposition in the last two layers of the HCAL. 
In this instance any tracks pointing to HCAL azimuthal boundaries were rejected. 

The identification criteria were studied with test samples of real data. The efficiencies 
of all the main selection cuts were tested using a sample of hadrons selected by tagging 
TT^'s in the HPC. This test sample allowed for an accurate calibration of all the main 
selection variables across the whole range of cos^* and cosip, the two variables used in 
the fits to the Michel parameters and the anomalous tensor coupling. 
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Remaining background from e^e^ -^ e+e^ and e^e~ — *> ^^fi^ events was suppressed 
by demanding that the particle in the opposite hemisphere to the identified hadron had 
a measured momentum of less than O.Spbeam- The four-fermion events e"'"e~ -^ e"'"e~r+r^ 
required no further suppression. 

A total of ~ 56000 r — >■ h{nT[^)vr candidates were selected from the data. The effi- 
ciency of selection within the 47r angular acceptance was 37%, the background arising from 
r 7^ h(n'K^)Pr processes was estimated to be (2.43 ± 0.73)% from e"'"e~ — >• e"'"e~ events, 
(0.40 ± 0.12)% from e+e" -^ /i+zi" events (0.10 ± 0.03)% and from e+e" -^ e+e-r+r" 
events (0.23 ± 0.07)%. 

6.5 The two-dimensional selection 

As described in Section ^, in order to measure the Michel parameters most efficiently it 
is necessary to use two-dimensional spectra. It was required that the events satisfied the 
preselection cuts and that there was one identified candidate r decay in each hemisphere. 
This therefore produces 20 (15 two-dimensional and 5 one-dimensional) distributions 
consisting of e/x, ee, /i/z, e/il,[] e/i2, e/i3, /i/il, /i/i2, /i/i3, /il/il, h2h2, hShS, hlh2, hlhS, 
h2h3, eX, fiX, hlX, h2X and h3X, where the two identified particles in each correspond 
to the two hemispheres in the event. The X in the event is an unidentified r decay with 
either one or three charged particles. In this case only the hemisphere with the identified 
track is used. 

In most of these channels it is required that the r preselection cuts be satisfied in 
order that non-r backgrounds be suppressed. This is not true for the e/i channel in 
which no preselection cuts were necessary as the external background required no further 
suppression. To suppress remaining cosmic ray background in the /i/x and the fiX samples 
it was required, in one-versus-one charged particle topologies, that at least one of the 
charged particle tracks extrapolated to within 0.3 cm in the r — plane of the interaction 
region. For the one-dimensional distributions, eX, fiX, hlX, h2X and h3X, the cuts 
to remove external backgrounds follow those already outlined in the previous sections 
describing the one-dimensional selections. 

The number of events selected, the efficiency of selection within the fiducial volume 
and momentum acceptance and the backgrounds can be seen in Tables |^ and ^ 



channel 


efficiency(%) 


z-J- 


-^T+T- 


-^ {evT'){^vV) 


72.95 ± 0.23 


Z"- 


^T+T- 


-^ {evT'){evV) 


50.43 ± 0.36 


Z^- 


^T+T- 


—>■ (/iZ/I7) (/iZ/I/) 


82.77 ±0.27 


z^- 


^T+T- 


-^ (ez/I7)(h(n7r°)z/) 


47.08 ±0.15 


z^- 


^T+T- 


-^ (/iz/I7)(/i(n7r°)z/) 


60.23 ±0.15 


z^- 


^T+T- 


-^ (/i(n7rO)z/)(/i(n7rO)z/) 


37.62 ±0.13 



Table 2: The efficiencies of selection in the angular and momentum acceptance for the 
two-dimensional analysis in the 1994 data set. The efficiencies were similar for the other 
years. The errors are purely statistical. 



^ where hl,h2 and h3 are hadrons in the invariant mass bins Mi„v < 0.3 GeV/c^, 0.3 GeV/c^ < Mj^v < 0.95 GeV/c-^ 
and Mi„v > 0.95 GeV/c^ respectively 
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T^T^ decay 


no. of candidate 


internal 


non-r+r 


modes 


events 


background(%) 


background(%) 


e-e 


1405 


6.01 ±1.80 


7.19 ±2.16 


e-/i 


3495 


4.43 ±1.33 


0.60 ±0.18 


e-hl 


1804 


5.31 ±1.59 


0.60 ±0.18 


e-h2 


3324 


3.92 ±1.18 


0.11 ±0.03 


e-h3 


1088 


4.01 ±1.20 


0.21 ±0.06 


e-X 


6377 


2.96 ±0.89 


4.66 ±1.40 


li-H 


2116 


2.61 ±0.78 


3.89 ±1.17 


fx-hl 


2160 


3.88 ±1.16 


0.30 ±0.09 


fi-h2 


4454 


2.13 ±0.64 


0.54 ±0.16 


fi-h3 


1480 


2.09 ±0.63 


0.80 ±0.24 


^-X 


8632 


1.58 ±0.47 


1.36 ±0.41 


hl-hl 


571 


4.89 ±1.47 


0.20 ±0.06 


hl-h2 


2271 


3.19 ±0.96 


0.10 ±0.03 


hl-h3 


730 


3.12 ±0.94 


0.15 ±0.05 


hl-X 


5104 


2.54 ±0.76 


1.87 ±0.56 


h2-h2 


2295 


1.69 ±0.51 


0.12 ±0.04 


h2-h3 


784 


1.80 ±0.54 


0.01 ±0.01 


h2-X 


9342 


0.93 ±0.28 


0.33 ±0.10 


h3-h3 


278 


1.93 ±0.58 


0.01 ±0.01 


h3-X 


3058 


0.94 ±0.28 


0.16 ±0.05 



Table 3: The number of selected events (column 2) and backgrounds (columns 3 and 4) 
for the selection described in the text. The backgrounds are quoted for the 1994 data set 
only. They were similar for the other years. A total of 60768 events were selected in the 
1992-1995 sample. 
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7 The extraction of the Michel parameters 

The values of the Michel parameters, p, 77, ^ and C,S together with the tau polarisation, 
Vt, and the tau neutrino helicity, h,^^, are extracted from the data using a binned maxi- 
mum likelihood fit to all the combinations of r — > eFgi/^, r -^ iJ'T'^i'r and r — > /i(n7r°)z/^. 
In splitting the hadron sample into three invariant mass bins one is left with 15 two- 
dimensional and five one- dimensional distributions where only one r decay has been 
exclusively identified in an event. 

The likelihood function is defined as: 

^=nn 'v ., (33) 

where n^- is the number of observed events in selected class c in the bin denoted by the 
indices i,j. The predicted number of events in this bin is af- and is given by 

«^ = ^^ E Vf ^r T^h + ^ + ^T^~^^- (34) 

The detector resolution matrix 7^^/^ gives the fraction of reconstructed signal events with 
generated fit variable in bin k' which are reconstructed in bin k. £?• describes the t^t^ 
selection efficiency as a function of the reconstructed fit variables in the two r decay 
hemispheres. The H. and E matrices were obtained from the full detector simulation. 
The matrix T contains the two-dimensional distribution corresponding to Eqn. |19|and the 
dependence on the fitted parameters. The construction of T, taking into account mass, 
radiation, and hadronic modelling effects, is described below. The number of background 
TT events per bin is 6^j^, and was not varied as a function of the fitted parameters. The 
non-rr background per bin, 6^!°°'^"^^, was normalised to the luminosity of the data. The 
signal and r background were then normalised keeping their ratio constant so that the 
integrals of the predicted fit distributions were the same as the total number of events 
seen in the data. 

This method accounted for correlations between the r"*" and r~ in an event arising 
from geometric detector reconstruction effects, described by the detector simulation, and 
physical effects such as longitudinal spin correlations, electroweak and QED corrections, 
described by the KORALZ program. Near the Z pole, photonic radiative effects are 
a strong function of the centre-of-mass energy. In the derivation of the efficiency and 
resolution matrices, t^t~ simulation samples have been used for the different centre-of- 
mass energies with proportions corresponding to the data sample. 

The hi[x) polynomials describing the leptonic decay spectrum shown in Fig. ^ do not 
take into account mass effects or radiative corrections. These effects were introduced 
by Monte Carlo methods using KORALZ and a modified version [^ of the TAUOLA 
program to generate distributions corresponding to the hi[x) polynomials. The TAUOLA 
program models leptonic r decays with the matrix element containing exact 0[oi) QED 
corrections. The modified version contained a generalisation of the Born level part of 
the matrix element which permitted the setting of non-Standard Model values for the 
Michel parameters. The part of the matrix element describing the QED corrections was 
calculated assuming V — A couplings. The part of the matrix element proportional to 
OiQED is small and it was assumed that for observed variations of the Michel parameters 
the change in the spectra due to changes in the radiative corrections could be neglected. 

For semi-leptonic r decays, the distributions were obtained from linear combinations 
of distributions generated with h^^ = — 1 and either positive or negative helicity states of 
the decaying r. 
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In the fit it was assumed that hy^ had the same value for all the semi-leptonic decay 
modes. 

In the fit assuming lepton universality the value of rj can be constrained using the mea- 
sured values of the leptonic branching ratios in Eqn. |Tl]. The branching ratio results pO[] 



were obtained from the DELPHI data in the years 1991 through 1995. The value of t] 
was constrained with the addition of the following quantity to the log-likelihood function 



^const _ l{r]-r]Br 



where rjBr is the value obtained from the leptonic branching ratio measurement and ArjBr 
is the error on this measurement. 

It must however be noted that obtaining Eqn. |Tl| involves an integration over the final 
state momenta, the implications of which have to be accounted for when setting a limit 
on 7] based on experimentally measured branching fractions. Since t] affects the shape of 
the muon momentum spectrum as well as the total decay rate, it is necessary to study the 
effect of the cutoff on the muon momentum identification which is at x'^ = p^/pheam = 0.05. 
As a function of the normalised laboratory muon momentum x = p/pbeam the number of 
events observed between momentum x and x + dx can be written as 

dN = No [a{x) + Kr]h{x)] dx. (36) 

By analogy with Eqn. H, the polynomial a{x) = Hq^x) + ^hp{x) — Vrih^ix) + ^h^s{x)] is 
the appropriate linear combination of polynomials for Standard Model couplings at LEP 
energies, while b{x) = hj^{x). The constants Nq and K can always be chosen such that 
the integrals of a{x) and h{x) over the whole momentum range are normalised to 1. If r/ 
is non-zero, the number of events observed would be 



Nobs = No 



nrfTnax rpffiax 

/ a{x)dx + Kt] / b{x)dx 



(37) 



The event generator used to compute the acceptance corrections, KORALZ/TAUOLA, 
assumes that t] equals zero. In other words, the branching ratio is derived assuming that 
the total number of r — >• /iz/^z/,- decays produced can be estimated as 

A^o'^* = Nobs X J (38) 

j^c a[x)dx 

where the integral is obtained from simulation. Hence, instead of correcting to obtain 
Nq^^ = Nq + Kt], the estimate of the corrected number of events becomes 



AT-* = No 



max 

r^ Lc b{x)dx 



J^c a[x)dx 



(39) 



The ratio between the integrals is readily calculated numerically by generating the 
full distribution in the r rest frame and boosting the momentum to the lab frame. It is 
found that the ratio between the integrals equals 0.96 when integrating from x*^ = 0.05. 
Ignoring effects due to r^ in r ^ eT/gZ/^ decays, the relation 



Br{T -^ ez^T-Fg) 



/rn,,\ m,, /rn,,\ 

f(^)+3M^g(^h,. (40) 



should be used to extract r]^ from the DELPHI tau leptonic branching ratios instead of 
Eqn. O. 
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Using the techniques outhned above together with background distributions obtained 
from the simulated data a six parameter and a nine parameter fit were performed, with 
and without the assumption of lepton universahty respectively, over a sample of ~60000 
r pair candidates. The one- dimensional projected distributions for each r decay class 
are shown in Figs. and ^, together with the fitted distributions obtained from the six 
parameter fit. 

A number of cross-checks were performed to check the stability of the result with 
respect to the selection cuts and binning effects. Each of the r^r^ preselection cuts was 
varied by 10% of its value (in the case of the Erad and Prad by 5% corresponding more 
closely to their resolution); no variation in the results was observed beyond those expected 
from statistical fluctuations. A similar process was performed for the main cut criteria 
in the selection of the different r decay modes classes. Again no variation was observed 
in the fit results beyond that expected from statistical fluctuations. The binning used to 
define the hadronic decay classes hi, h2 and h3 was varied; no unexpected variation in 
the results was observed. 

The systematic effects studied are decribed below and summarised in Tables ^ and ^. 

One source of systematic uncertainty arose from the finite amount of simulated data 
available. 

An uncertainty due to the r branching ratios was obtained by varying the branching 
ratios by the uncertainties on the world average values, repeating the fit and taking 
the change in the result as an estimate of the systematic uncertainty. Conservatively, 
the background from r decays and the t~^t~ backgrounds were varied by 30% and the 
change in the results of the fit taken as an estimate of the systematic uncertainty. 

The dependence of the selection efficiency on the fit variables for the different r decay 
modes was studied using data test samples or by redundancy of the different detector 
components as described in each of the relevant sections. The resulting systematic un- 
certainties were estimated by varying the selection efficiency in the simulation as a linear 
function of the fit variable. The magnitude of the variation was taken from the statistical 
uncertainty on the gradient derived in a straight line fit to the ratio of the measured 
efficiencies in data and simulation as a function of the fit variable. 

Systematic uncertainties were attributed for detector calibration effects. The charged 
particle reconstruction momentum scale was varied by its uncertainty, the analysis re- 
peated, and the resultant variation in the results taken as an uncertainty. This in particu- 
lar affected the muon channel parameters. The effects due to knowledge of the momentum 
resolution were also taken into account but were much smaller. 

The electron momentum estimator pei was calibrated on data using both radiative and 
non-radiative Bhabha events. Its scale was calibrated with a precision of 0.5%, limited 
by the statistics of the data test samples. The systematic uncertainties on the various 
fit parameters were estimated in the same way as for the momentum scale. A smaller 
contribution arose from the knowledge of the resolution on p^i- 

The neutral electromagnetic energy scale was known with a precision of 0.2%. The 
related uncertainties were estimated in an analogous manner to those due to the momen- 
tum and Pel scales. This affected mostly the hadronic decay modes with vr^'s in the final 
state. 

The uncertainty in the energy scale in the HCAL had a negligible effect, as did the 
uncertainty in the efficiency of the muon chambers. 

The systematic uncertainty contribution arising from the HCAL response to hadronic 
showers was estimated by varying within their statistical errors the corrections to the 
shower penetration in the simulation taken from the data test samples. 



22 



The hadronic invariant mass scale uncertainty is dominated by the neutral energy and 
charged particle momentum scale. Any serious discrepancy between simulation and data 
would be evident in the hadronic invariant mass distribution such as that shown in Fig ^, 
where agreement is good. Additional checks have been made on the spatial resolution of 
the electromagnetic showers in the HPC These effects were found to be small compared 
with those due to energy scale and resolution. 

Imperfections in the modelling of the photon reconstruction efficiency could lead to a 
poorly modelled cross-talk between the different invariant mass classes in the inclusive 
hadronic selection as well as affecting the reconstruction of the 9* and ip angles used 
in the fit. From a study |3ll] of various distributions related to reconstructed photons, 
such as multiplicity and energy distributions, it was estimated that the neutral photon 
reconstruction efficiency was known to better than 4% averaging over the whole of the 
HPC taking into account dead space and threshold effects. The systematic uncertainty 
attributed to this was estimated by randomly rejecting 4% of photons and the change in 
the results was included under the heading calibration. Further cross-checks of the HPC 
reconstruction in the inclusive hadronic sample included reclassifying energy depositions 
associated to the charged particle track as neutral particles. This had a negligible effect 
on the results, indicating that both the mis-association of photon and 7r° showers to the 
charged hadron and the description of the hadronic interactions associated to the hadron 
charged track and misidentified as electromagnetic showers were well described by the 
simulation. 

The uncertainties due to radiative corrections in hadronic r decays and modelling of 
the ai have been estimated to give a systematic uncertainty of 0.001 on Vr for the inclusive 
hadronic polarisation analysis in [^. This has been included as a systematic on Vr and 
the systematic uncertainty has been propagated through to the other fit parameters. 

Fig. shows discrepancies between data and the fitted distributions for the r -^ hlu^. 
sample. Studies of the quantities used to select the sample exhibited no obvious effect 
which could account for this. A cross-check, performing the fit excluding the data from 
the discrepant regions (the ranges [-0.9,-0.7] and [0.8,1.0]) of the r — > /ilz/^, showed 
variations which were consistent with statistical fluctuations. Conservatively, systematic 
uncertainties were estimated for this effect by taking half of the variation in the fit 
results when forcing the fits to go through the data points in the quoted ranges. These 
are included as a contribution to the "calibration" uncertainty in Tables |^ and ^. 
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Vr 
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i5 


Kr 


MC stats 
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0.0035 
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0.0104 


0.0103 


0.0039 


rBR's 


0.0002 


0.0006 


0.0014 


0.0004 


0.0012 


0.0020 


Backgrounds 


0.0251 


0.0115 


0.0011 


0.0030 


0.0126 


0.0093 


Efficiency 


0.0005 


0.0023 


0.0037 


0.0013 


0.0027 


0.0014 


Calibration 


0.0144 


0.0146 


0.0065 


0.0281 


0.0229 


0.0034 


Decay modelling 


- 


0.0009 


0.0010 


0.0007 


0.0004 


0.0010 


7] const. 


0.0232 


0.0070 


- 


- 


- 


- 


Total Syst. 


0.037 


0.020 


0.008 


0.030 


0.028 


0.011 


Statistical 


0.036 


0.023 


0.012 


0.070 


0.070 


0.027 



Table 4: The systematics on the parameters for the six parameter fit with the assumption 
of universality. The statistical error is shown for comparison. 
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Figure 7: The projections of the fitted distributions for the six parameter fit for the two 
fully leptonic decay channels and the semi-leptonic candidates from the lowest invariant 
mass bin. The line is the result of the fit, the points are the data and the shaded area is 
the sum of the backgrounds. 



The six parameter fit assuming lepton universality (including the constraint on t] from 
the leptonic branching ratios) gave the following results: 



V = 


= -0.005 ±0.036 ±0.037 


p = 


= 0.775 ±0.023 ±0.020, 


e = 


= 0.929 ±0.070 ±0.030, 


^s - 


= 0.779 ± 0.070 ± 0.028, 


^Ur = 


= -0.997 ±0.027 ±0.011 


Pr = 


= -0.130 ±0.012 ±0.008 



The results were found to be stable as a function of the year of data taking. After cor- 
recting for effects including the photon propagator and the -y/i dependence, the result on 
Vt obtained in this analysis can be compared to the polarisation parameter Ar obtained 



in the dedicated analysis of the DELPHI data [^. Correcting Vt as described in |^l]] 
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Figure 8: The projections of the fitted distributions of the w variable described in |[T5| 
from the six parameter fit for the semi-leptonic candidates selected in the second and 
third invariant mass bins. The line is the result of the fit, the points are the data and 
the shaded area is the sum of the backgrounds. 



the result At = 0.134 ±0.014 is obtained for this analysis. This is in excellent agreement 
with the result A = 0.1359 ± 0.0096 from |3l|. 

The x^ /Ndof for the fit was 1984/2009. The parameters are correlated and the corre- 
lation matrix is given in Table |^. 

The variable P^, defined in Eqn. |12|, represents the probability of a right-handed r 
decaying into a lepton of either handedness. This was calculated to be 



Pr 



-0.038 ±0.066 ±0.029. 



A one-dimensional ffi to t] was also performed. In setting the other Michel parameters to 
their Standard Model values and applying the branching ratio constraint (Eqn. |40| ) the 
value of 7] was found to be 



7] 



-0.009 ±0.033 ±0.024. 





^M 


Pe 


Pm 


Vr 


^e 


i. 


^e^e 


i^.5^ 


Kr 


MC stats 

rBR's 

Backgrounds 

Efficiency 

Calibration 

Decay modelling 


0.047 
0.003 
0.138 
0.010 
0.039 
0.002 


0.0054 
0.0008 
0.0230 
0.0034 
0.0278 
0.0010 


0.0144 
0.0016 
0.0414 
0.0017 
0.0076 
0.0004 


0.0018 
0.0016 
0.0013 
0.0036 
0.0069 
0.0010 


0.0177 
0.0017 
0.0058 
0.0020 
0.0438 


0.028 
0.001 
0.044 
0.003 
0.018 
0.001 


0.019 
0.001 
0.018 
0.001 
0.034 
0.001 


0.019 
0.001 
0.018 
0.001 
0.033 
0.001 


0.0039 
0.0003 
0.0093 
0.0014 
0.0045 
0.0010 


Total Syst. 


0.15 


0.037 


0.045 


0.008 


0.05 


0.06 


0.04 


0.04 


0.011 


Statistical 


0.32 


0.036 


0.098 


0.012 


0.12 


0.19 


0.12 


0.13 


0.028 



Table 5: The systematics on the parameters for the nine parameter ffi without the 
assumption of universality. The statistical error is shown for comparison. 
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Table 6: The correlation matrix for the six parameter fit. 



The nine parameter fit without any assumption of lepton universahty gave the follow- 
ing results: 



Pe 
P. 



0.72 ±0.32 ±0.15, 
0.744 ±0.036 ±0.037, 
0.999 ±0.098 ±0.045, 
1.01 ±0.12 ±0.05, 
1.16 ±0.19 ±0.06, 
0.85 ±0.12 ±0.04, 
0.86 ±0.13 ±0.04, 
-0.991 ±0.028 ±0.011, 
-0.131 ±0.012 ±0.008. 



The parameters are correlated and the correlation matrix is given in Table |^. 
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^e^e 
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0.087 


^f^S^ 
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Table 7: The correlation matrix for the nine parameter fit. 



The values of the Michel parameters for the process r -^ ^iV^Vr are less precisely known 
than those from the r — > eV^^Vr channel. This is because for the r -^ ^iV^Vr channel one is 
also measuring the r] parameter which has all its sensitivity in this channel. The 77^ and 
p^ parameters are both at the level of ~ 2a away from the Standard Model predictions. 
These two parameters are very highly correlated. In setting rj^ to its Standard Model 
prediction value of one obtains the following results: 

?7^ = (fixed), 

Pe = 0.755 ±0.036 ±0.037, 

p^ = 0.789 ±0.028 ±0.012, 
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1.00 ±0.12 ±0.05, 
0.87 ±0.11 ±0.03, 
0.86 ±0.12 ±0.04, 
0.733 ±0.094 ±0.030, 
-0.995 ±0.028 ±0.011, 
-0.129 ±0.012 ±0.008. 



The presented measurements show no deviations from the predictions of pure V — A 
couphngs in r decays. 

As mentioned in Section |^ the Michel parameters are restricted by boundary condi- 
tions. The physically allowed regions for various pairs of the parameters p, ^ and ^6 are 
shown in Fig. ^ along with the experimentally determined values for these parameters, 
for the fit assuming lepton universality. In forming the contours the likelihood function 
is minimised with respect to the other four parameters in the fit. One can see that the 
contours enter into the disallowed regions due to the finite experimental resolution. The 
disallowed regions in Fig. |^ are in fact dependent on three of the Michel parameters. The 
disallowed regions shown are presented with the Michel parameters set at their Standard 
Model values for simplicity. These regions will therefore move around to encompass more 
of the fitted contours if the Michel parameters are set at their measured values. 
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Figure 9: The contours corresponding to {logL ± y), where i = 1,2, 3, for the six param- 
eter fit. In forming the contours the likelihood function is minimised with respect to the 
other four parameters in the fit. In each plot the allowed region corresponds to the case 
where all the Michel parameters are fixed to their Standard Model values. 
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8 Extraction of the coupling n^ 

The spectra of the r decay products were used to extract the parameter k^ . To 
estimate the theory prediction of the spectra distortion in the case of k^ ^ the Stan- 
dard Model simulated data were used with the events re-weighted in the following way. 
For the generated values of the r helicity and the final lepton momentum, the value 
dT/dxi{xi, K^) was calculated according to Eqns. ^ and ^ The ratio dr/^x(x'^o) ^^^ 
then used as an event weight to produce the simulated spectrum with non-zero tensor 
coupling. In the case of the r multipionic decays Eqn. |26| was used to generate event 
weights. The radiative corrections to r production and r decay were taken into account 
by the KORALZ 4.0 and TAUOLA 2.5 programs and the variable xi was defined using 
the lepton energy after all radiation. It was assumed the effect of the tensor coupling was 
small and that the effects of radiative corrections to the tensor coupling contributions 
could be neglected. 

The value of the tensor coupling parameter was then extracted from a log likelihood 
fit of the simulated spectra to the real data, with k}^ as a fit parameter. One-dimensional 
spectra of xi were used in the case of leptonic r decays and the two-dimensional spectra 
of (cos^*, cos^^) for semi-leptonic decays. To increase the sensitivity of the semi-leptonic 
channel further, the region of reconstructed invariant mass between 0.3 and 1.7 GeV/c^ 
was divided into five bins and the fit was performed in each bin simultaneously. This 
reduced the statistical error of the fit by about 10%. The region of invariant mass below 
0.3 GeV/c^ was not used because it was dominated by r ^ nur decays which have no 
sensitivity to the tensor coupling. 



The illustration for the channel r -^ ^j^^i't is given in Fig. |T0| which shows the dif- 
ference between the real data and the Standard Model simulated data prediction. Also 
shown is the difference between the best fit simulated data and the Standard Model sim- 
ulated data. The sample used for this channel was an equivalent data set to that used for 
the Michel parameter analysis together with r — * fiV^Ur decays detected in the end-cap 
region of the detector. The selection follows that described in [^ . 

The systematic uncertainty for the muon channel received contributions from the lim- 
ited simulation statistics, the dimuon background level and the calibration of the mo- 
mentum scale of the charged particle track reconstruction. For the eUeUj. final state the 
main contributions to the systematic uncertainty are the limited simulation statistics, the 
level of the Bhabha background and the calibration of the electron energy estimator pei. 
For the hadronic selection the main systematic is the calibration of the neutral electro- 
magnetic shower energy scale. Other significant contributions arise from the efficiency of 
photon detection in the HPC and the knowledge of the resolution on the fit quantities 
6* and ip. The calibration of the momentum scale for the charged hadron also gives a 
contribution. Other sources of systematic uncertainty common to some or all channels 
are the modelling of the momentum resolution and the electromagnetic energy resolution, 
and the uncertainty on the r branching ratios. 

The results of the fits for different decay channels were the following: 

K^ = +0.162 ± 0.078 ± 0.030, 



r — > ez/eZ/r 

,0\ 



-W 



= -0.043 ±0.057 ±0.032, 
K^ = -0.122 ± 0.059 ± 0.025. 



r -^ h{mi )i>r 
Combining these taking into account correlations gave 



n^ = -0.029 ± 0.036 ± 0.018, 
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Figure 10: An illustration of the fit to the tensor coupling parameter using the decay 
r -^ jjlV^i't-. Upper plot: spectra of the normalised muon momentum for data (points 
with error bars), background (black) and the best fit simulated data (solid line). Lower 
plot: the difference between the measured spectrum and the Standard Model prediction 
(points with error bars); the solid line shows the difference between the best fit simulation 
and the Standard Model simulation. 



where the first uncertainty is statistical and the second is systematic. The x^/^rfo/ for 
this combination is 7.6/2 corresponding to a probability of having a worse x^ /N^of of 
2.3%. The measured value of the anomalous tensor couphng, n^ corresponds to a 90% 
allowed interval of -0.096 < k^ < 0.037. 



9 Interpretation of the results 

The results of this analysis can be used and interpreted in a number of different ways. 
It is interesting now to use the measured values of the Michel parameters to explore these 
avenues and investigate the possible existence of new physics beyond the Standard Model. 

The measured values of the parameters ^ and ^5 were used to estimate the probability 
P^ of a right handed tau decaying into a lepton of either handedness. Following the 
technique outlined in ||32|, and only allowing the value of P^ to be between and 1, the 
corresponding upper limit on this quantity was found to be: 



P^ < 0.081 



at 90% C.L. 



Using the constraints described in Section ^ one can set limits on the coupling constants 
by forming positive definite expressions from the measured parameters. The best limits 
on gf^^ and gj^j^ are derived from Eqn. |12] using the limit on P^. Eqn. |16] constrains 
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g^j^, while the best constraints on (7^^ and g^ifi are derived from Eqn. |T^. The best 
hmits on the f^^^^, gf^j^ and (7^^ are obtained from Eqn. 13. Absorbing the freedom given 



by the unknown overall phase, the coupling gY^ is taken to be real and positive. Only 
the coupling gl^^ cannot be constrained as one cannot distinguish between ^ff^ and the 
Standard Model coupling g^L- It could be constrained by measuring the cross-section for 
inverse r decay ^ . The overall normalisation factor A in Eqns. |^ and ^ was fixed to a 
value of 16. 

The 90% confidence level upper limits on the coupling constants, derived as described 
in 1^, are given in Table H for the fits with and without lepton universality. The 
parameters gf^^^ and g^^^ are coupled together in Eqn. 0. The limits obtained for (7^^ 
assume gf^^ = 0; if this condition is relaxed, the limits obtained for g^j^ are poorer 
than the normalisation constraint. The results of the fit assuming universality are also 



illustrated pictorially in Fig. 11 
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Table 8: 90% CL. upper limits on the magnitudes of the complex coupling constants. 
The 2nd column contains the results assuming e-/i universality. The 3rd and 4th columns 
display the results for the electronic and muonic decay modes respectively. The fifth 
column shows the maximum physically allowed value for the parameter. 



In its minimal version, the Higgs mechanism is implemented by adding only one dou- 
blet of complex scalar fields resulting in one additional physical scalar state, electrically 
neutral, commonly referred to as the standard Higgs boson. One can postulate extensions 
to this by adding, for example, one more doublet of complex scalar fields, which leads 
automatically to five physical states (three neutral and a pair of charged Higgs bosons), 
after the spontaneous breaking of the SU(2)l x U(1)y symmetry to give mass to the W^ 
and Z gauge bosons. These charged Higgs bosons contribute to the r — > Wi^r, l=e,fi 
decay through a scalar coupling, given at Born level by 



9i 



m/mT-tan^/3 

2 ' 



(41) 



for negligible neutrino masses |1T0| , |33| , |34| . Tan/5 is the ratio of the vacuum expectation 
values of the two Higgs doublets and mH± is the mass of the charged Higgs boson. For 
left-handed neutrinos the couplings are of the type gfui, and the Michel parameters rji 
and ^i, (with l=e,fj,), can be written in terms of gfifif 



m 



9'rr,i/2 
- {9'rrM 



(42) 
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Figure 11: Limits (at 90% CL) on the coupling constants from the Michel parameter fit 
assuming lepton universality. The dark shaded areas are the allowed regions. 



1 - fe.^/2)^ 
1 + (^L,//2)^ 



(43) 



while 



4' 



With these relations mii± can be extracted using the measured values 



of the Michel parameters. Using the results presented in this paper one obtains a lower 
limit on the mass of the charged Higgs of 



mH± > 1.17 X tan/? GeV/c^ at 90% C.L. 



This limit is not competitive with those from direct searches, unless tan/? has an unex- 
pectedly large value. 

Another extension to the Standard Model which can be related to the Michel pa- 
rameters involves the postulate that parity violation is caused by spontaneous symmetry 
breaking. A question which arises in the Standard Model is why the doublets are left 
handed and the singlets right-handed. Left-right symmetry implies that the Lagrangian 
is both charge and parity invariant before spontaneous symmetry breaking and that CP 
violation arises due to the non-invariance of the vacuum. These left-right symmetric mod- 
els assume the existence of a second pair of W bosons, and the weak eigenstates W^l,/? 
are mixtures of the mass eigenstates W\^2 [p5| , p6| . One can introduce the mass ratio^ a 
of the mass eigenstates, 

a = m?,,/mL, (44) 



Wxi 



the limits presented for the two doublet and left-right symmetric models correspond to the value at which the likelihood 
has dropped to give the corresponding measure of confidence. 

^this mass ratio is commonly referred to as /3 in the literature. In order to avoid confusion with the parameter introduced 
above, tan/3, the name has been changed to a 
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and the mixing angle ( between the weak and mass eigenstates. The Michel parameters 
p and ^ take the form 

3 A J^ . . A . . _4a 

r 

l-a2 



p = ^cos^Cf^l + tan^C+Y^^tan^C), (45) 



e = cos^C(l-tan^C)T— ^, (46) 

while rj = and 5 = |. The z/^ polarisation parameter in r hadronic decays takes the 
form IB? 



K^ = -1 + 2C' + a(4C - 8C') + a'(2 - 8C + 12C'). (47) 

Using these relations, the measured value of mwi and the measured values of the Michel 
parameters, one can place limits on 171^2 and (, with the caveat that the right-handed 
neutrino must be light enough compared with the r to be produced without kinematical 
suppression. Taking the measured values of the Michel parameters and mwi = (80.41 ± 



0.10) GeV/c^ |3^ gives the following limits on rriWi and the mixing angle (: 

mw2 > 189 GeV/c^ at 90% C.L.; 
-0.141 < C < 0.125 rad at 90% C.L. 

The 68%, 95% and 99% confidence level contours on the C-^tT'W2 plane are shown in 
Fig. |T2|. The x^ function exhibits a slight minimum at mvK2 = 290 GeV/c^, ( = —0.01. 



The distribution then exhibits an allowed region extending to infinite mw^ where the 
change in x^, for (^ = 0, is 0.13 compared with the minimum. For the case of no mixing 
(^ = 0, and the lower limit on the W2 mass becomes 

mw2 > 204 GeV/c^ at 90% C.L.. 



10 Conclusions 

A precise measurement of the Michel parameters and the z/,- helicity has been presented, 
together with limits on the anomalous tensor coupling. 

A simultaneous fit to the Michel parameters and the Uj. helicity assuming e-/i univer- 



sality and using the DELPHI exclusive leptonic branching ratio measurements [30| gave 
the following results: 

rj = -0.005 ±0.036 ±0.037, 

p = 0.775 ±0.023 ±0.020, 

^ = 0.929 ±0.070 ±0.030, 

^S = 0.779 ± 0.070 ± 0.028, 

h^^ = -0.997 ±0.027 ±0.011. 

A fit to the Michel parameters and the u^- helicity not assuming universality gave the 
following results: 

rif, = 0.72 ±0.32 ±0.15, 
Pe = 0.744 ±0.036 ±0.037, 
p^ = 0.999 ±0.098 ±0.045, 
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Figure 12: Contours corresponding to 68%, 95% and 99% confidence levels on the ( 
versus mw2 plane for the fit to the left-right symmetric model. 



^e = 1.01 ±0.12 ±0.05, 

^^ = 1.16 ±0.19 ±0.06, 

^Je = 0.85 ±0.12 ±0.04, 

^i^Sf, = 0.86 ±0.13 ±0.04, 

h^^ = -0.991 ±0.028 ±0.011. 

In both fits the average r polarisation was left as a free parameter. 

The world averages for the z/^ helicity and for the Michel parameters p, ^ and ^5, 
both with and without the assumption of lepton universality, are dominated by the re- 



sults ||3^ from the CLEO experiment. The results presented here are consistent with the 
world averages and have a precision between 1.5 and 2.5 times poorer than the CLEO 
measurements for these parameters but are of a similar or higher precision than other 
measurements, taken at LEP [|,|0|,|T],|2l or by ARGUS g. 

Fixing h^^ to the standard model value of —1 and constraining Vt from other mea- 
surements of the effective weak mixing angle would give a reduction in the quoted errors 
on the Michel parameters due to the correlations present in the simultaneous fits to all 
the parameters. 

The measurement of r] presented here is the most precise recorded to date, due to the 
combination of the precise leptonic branching ratio measurements and the measurements 
of the spectra. 

A measurement of the tensor coupling k^ has been performed for the first time in r 
decays, yielding the result 



K, 



.w 



-0.029 ±0.036 ±0.018, 



consistent with zero. 
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The presented results are consistent with the Standard Model and limits have been 
placed on the magnitudes of the complex coupling constants. The V — A assumption 
is however still not fully verified. Future results from B factories, complemented by a 
measurement of inverse r decay, will allow a full determination of the Lorentz structure 
of the r. 
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